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Abstract 

We consider nuclear spin-lattice relaxation rate resulted from a diffusion equation 
for rotational wobbling in a cone. We show that the widespread point of view that 
there are no analytical expressions for correlation functions for wobbling in a cone 
model is invalid and prove that nuclear spin-lattice relaxation in this model is exactly 
tractable and amenable to full analytical description. The mechanism of relaxation is 
assumed to be due to dipole-dipole interaction of nuclear spins and is treated within 
the framework of the standard Bloemberger, Purcell, Pound - Solomon scheme. 
We consider the general case of arbitrary orientation of the cone axis relative the 
magnetic field. The BPP-Solomon scheme is shown to remain valid for systems with 
the distribution of the cone axes depending only on the tilt relative the magnetic field 
but otherwise being isotropic. We consider the case of random isotropic orientation 
of cone axes relative the magnetic field taking place in powders. Also we consider 
the cases of their predominant orientation along or opposite the magnetic field and 
that of their predominant orientation transverse to the magnetic field which may be 
relevant for, e.g., liquid crystals. Besides we treat in details the model case of the 
cone axis directed along the magnetic field. The latter provides direct comparison 
of the limiting case of our formulas with the textbook formulas for free isotropic 
rotational diffusion. The dependence of the spin-lattice relaxation rate on the cone 
half-width yields results similar to those predicted by the model-free approach. 

Key words: NMR, spin-lattice relaxation, rotational diffusion, wobbling in a cone, 
model-free approach. 
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1 Introduction 



Diffusometry and relaxometry is a traditional and well developed branch of 

nmr CP, 0, [a, m, 0, 0, 0, is], in], ma, nu, na, usi, nu, usi, ng, 

HB], CI], US], |2Q], [21], [22], [23], I21], [IS], [26j. It is successfully applied to 
many complex systems such as proteins and homopolypeptides [2], [3], [27], 
U [29], [30], [15], [31], [32], [33], [31], [35], [36], [31], tissues [38], liquid 
crystals [9], [32], [15], liquids in porous glass [5J, [15], [IS] , polymers [15] . etc. 
Purely empirical way to interpret relaxation behavior in complex systems is 
to introduce a distribution of the correlation times, e.g., Cole-Davidson or 
Cole-Cole ones (see e.g., [8], [3D], [31], [32] and refs. therein). Usually this is 
done within the framework of the so-called model-free approach [33], [33], [E], 

[27J, [23], n M- 

At the same time many important theoretical issues still remain open. For 
instance free isotropic rotational diffusion is well investigated in details [17] . 
However it is applicable to a limited number of cases. In practice one is encoun- 
tered as a rule with some sort of restricted rotational diffusion, e.g., wobbling 
in a cone with half-width # treated rigorously by Wang and Pecora [33] • The 
latter model is widely used for interpreting NMR relaxation data pQ , [32] , [31] , 
[50], [51]. It has become a widespread point of view that there are no analytical 
expressions for correlation functions for wobbling in a cone model [50], [51], 
[39]. For instance in [5T] one can read: "although < P 2 (/t(0) • fi{t)) > cannot 
be evaluated analytically within the diffusion in the cone model, Analo- 
gous statement is reiterated in |49j:" Although there is no analytical expression 
for G m (t) in the case of this diffusion in a cone model, We do not agree 
with this point of view and the aim of the present paper is to show that nuclear 
spin-lattice relaxation in this model is exactly tractable and amenable to full 
analytic description. Here it is pertinent to recall a quotation from |51j:"To 
evaluate G m (t) exactly, one must solve time- dependent rotational diffusion 
equation subject to reflecting boundary conditions at 6$. Using this approach 
(Wang and Pecora [33]), one can express G m (t) as an infinite sum of exponen- 
tials, with amplitudes and time constants which are not closed-form functions 
of ^o"- Then the authors of [5T] circumvent the problem by constructing a 
simple but accurate closed- form approximation to G m {t). These results are 
sometimes referred to as quasi-exact [52] while the approximation is referred 
to as multi-exponential in contrast to the mono-exponential obtained in |50j . 
In the present paper we show that direct and stringent way to tackle the 
problem is nevertheless quite feasible. 

Thus in our opinion the quoted above assertion means not that statements 
like "there is no analytical expression for G m (t) in the case of the diffusion in 
a cone model" are valid but rather that exact analytic solution of the prob- 
lem merely has not been obtained yet. In this regard the notion of analyticity 
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should be refined. At the time the quotation from [51] was written only in- 
complete and erroneous Pal's tables for the values of v™ (the roots of the 
derivative of the associated Legendre functions with respect to the degree for 
the cone half- width O ) were available. This fact indeed made the expressions 
containing v™ to be "not closed-form functions of do", the formulas to be 
practically intractable and as a consequence the resulting theory to be non- 
analytic. However the situation was dramatically improved shortly after the 
publication of [51]. The appearance in 1986 of the Bauer's tables for the values 
of v™ as known functions of cos 9q [53] makes the problem of treating v™ con- 
ceptually equivalent to that of any habitual function. One should consider v™ 
as tabulated functions of 9q just as for instance a trigonometric function of 6q. 
In this sense a formula containing v™(Qq) since 1986 is as well analytic as that 
containing, e.g., sin(# )- In our opinion the issue about analytic tractability 
of the correlation functions within the framework of the diffusion in a cone 
model has long been in need for revaluation. 

In the present paper we precisely follow the instruction from the above quo- 
tation from [51]. We obtain analytic expressions for the exponentials from the 
infinite sum. Then with the help of a computer and Bauer's tables the problem 
of calculating the G m (t), the corresponding spectral densities and finally the 
spin-lattice relaxation rate becomes a routine though tedious procedure. We 
carry out these calculations for representative set of model parameters and 
various distributions of the cone axis relative the laboratory fixed frame (ex- 
ternal magnetic field) and plot spin-lattice relaxation rates. The dependencies 
obtained exhibit rich variety of interesting non-monotonic behavior. Our rig- 
orous quantum-mechanical treatment has a lucid classical analogy. The latter 
suggests some physical interpretation of the phenomenon. Also we find that 
our results are similar to those predicted by the model-free approach. Thus 
in our opinion the problem of calculating the correlation functions and spin- 
lattice relaxation rates within the framework of diffusion in a cone model gets 
full analytic solution. The same is true for spin-spin relaxation rate though 
this value is not addressed in the present paper to save room. 

The approach is discussed mainly for practically important case of an isolated 
two-spin system comprising a hetero-nuclear pair of non-identical spins, e.g., 
15 N — H in protein backbone or 13 C — H in protein side chains. However 
the results for homo-nuclear spin pair are also presented for the sake of com- 
pleteness. For wobbling in a cone we can extend the cone half-width 6q up to 
the limit of isotropic rotation 6q — > n to verify the coincidence of the results 
obtained with the known formulas. Only for the artificial potential of wob- 
bling in a cone the equation for restricted rotational diffusion can be solved 
in a stringent way without any approximation. For no other non-trivial model 
such exact treatment is possible. Thus the results of rigorous treatment of the 
rotational diffusion equation for wobbling in a cone can serve as a touchstone 
for approximations by necessity invoked to in the case of more realistic poten- 
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tials, such as, e.g., a harmonic one. In this regard the wobbling in a cone model 
is distinguished in the universe of models for restricted rotational motion. All 
other models require severe approximations and so long as the wobbling in a 
cone model is treated approximately it is merely one among many others. But 
as soon as this model is treated rigorously it takes an outstanding position 
of the exactly solvable one. Namely this fact motivates our attempt to revisit 
the investigation of the thirty years old problem of applying the Wang-Pecora 
model to NMR. 

We discuss the possible manifestation of the deviations of rigorous results 
for the spin-lattice relaxation rate from those of the approximate approach 
|43j . [5T] in the experiment. The requirements for the possibility of observing 
the above mentioned non-monotonic behavior in the experiment are found 
to be rather severe. Thus we anticipate that it can manifest itself provided 
it is sought purposefully and special concern is paid to the requirements. 
In these particular experimental situations the rigorous quantum-mechanical 
treatment developed in the present paper is expected to be more accurate com- 
pared with the approximate one of [13], [51]. By present no such experiments 
have been carried out. The rigorous treatment of the diffusion in a cone model 
enables one to reveal very subtle peculiarities in the behavior of the spin- 
lattice relaxation rate. That is why to verify them reliably in the experiment 
its special settings up should be deliberately devised. The advantage of the 
present approach compared with the approximate one lies mainly in the fact 
that the resulting theoretical description is completely coherent, i.e., results of 
the quantum-mechanical treatment agrees with those of its classical analogy 
and with those of [13], [51]. In our opinion such full picture enables one to 
gain more penetrating into the diffusion in a cone model and its possibilities 
to describe nuclear spin-lattice relaxation. 

The paper is organized as follows. In Sec. 2 the diffusion equation for restricted 
rotational wobbling in a cone is used to derive the joint probability density 
function. In Sec. 3 the latter is used to obtain the spectral densities of cor- 
relation functions for dipole-dipole interaction within the framework of the 
standard Bloemberger, Purcell, Pound (BPP) - Solomon scheme. Sec. 4-5 deal 
with the particular case that the cone axis is directed along the magnetic field. 
In Sec. 4 the spin-lattice relaxation rate for hetero-nuclear spin pair is obtained. 
In Sec. 5 that for homo-nuclear spin pair is obtained. In Sec. 6 the general case 
of arbitrary orientation of the cone axis relative the magnetic field is consid- 
ered. These results are applied to the case of isotropic random orientation 
(unweighted average) of cone axes relative the laboratory fixed frame. Also 
a model example of predominant orientation of cone axes along or opposite 
the magnetic field and that of their predominant orientation transverse to the 
magnetic field are considered. In Sec. 7 the results are discussed and the con- 
clusions are summarized. In Appendix A some known mathematical formulas 
are collected for the convenience of the reader. In Appendix B some technical 
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details of calculations are presented. In Appendix C the classical analogy of 
the quantum-mechanical treatment is considered. 



2 Rotational diffusion in a cone 



We choose a laboratory fixed frame so that its z axis of the Cartesian frame 
x, y, z is directed along the constant magnetic field. The random functions 
of relative positions of two spins specified below (see (THj) ) are defined in the 
corresponding spherical frame 9, <ft given by the polar angle 9 (counted from 
the z axis) and azimuthal one <j). We consider a general case that the cone axis 
is tilted at an angle ip relative the magnetic field. We direct the z' axis of the 
dashed Cartesian frame x', y', z' along the cone axis. The correlation function 
for wobbling in this cone we define in the corresponding spherical frame 9', 4>' 
given by the polar angle 9' (counted from the z' axis) and azimuthal one 
(f)'. Following [18] we consider a rod with the orientation specified by a unit 
vector u directed along its axis with spherical polar coordinates Q' = (9',<p'). 
In accordance with [47J we following Debye assume that the rotation of the 
rod can be considered as that of the hard sphere with radius a (a is the 
length of the rod) in a medium with viscosity rj. For the ordinary diffusion in 
a cone model the rod is allowed to diffuse freely within an empty cone with 
a maximum polar angle 9' = 9 . The symmetry axis of the cone is taken to 
be the z' axis. For the diffusion in a cone model, the polar angle is restricted 
(0 < 9' < 6> ) but the azimuthal angle is not (0 < ft < 2vr). 

Our aim is to consider ordinary diffusion for rotational motion in a cone. The 
probability density for finding the rod oriented in u at time t, i. e., ^(it,t), 
obeys the DE for rotational motion 

^M = D A*(M) (1) 



where D is the diffusion coefficient (DC) for rotation and A is the angular 
part of the Laplace operator in polar spherical coordinates 



1 

sin 2 9' 
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The DC for rotation D has the dimension cm? / s and is given by the Stokes 
formula 

87ia 3 r] ^ 
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where ks is the Boltzman constant and T is the temperature. 
Following [18] we write the solution of ([T]) as follows 

= 

oo oo 

E E exp[-CK m + l) J DUI]n? ) *(^(0))^ ) (^W) (4) 



n=l m=— oo 



where Y^m^ is the generalized spherical harmonics of degree v™ [48] . the 
symbol * indicates the complex conjugate and the values of v™ are determined 
by the boundary conditions on 9' defined by our diffusion in a cone model. 
The boundary condition says that there is no net change of the probability 
density at the boundary of the cone, i.e., 



89' 



(5) 



The values v™ are known functions of cos#o [H], [53]. They satisfy the re- 
quirement v~ m = v™ . The index n is defined such that v" 1 < v™ < v™ < ... 
for any m. The detailed calculations of v™ are presented in the tables [53] . 
The values of U™ for n > 1 increase sharply with the decrease of the confining 
volume. 

The solution (J1J is subjected to the initial condition 

0) = 5 (Of - fi'(0)) = 5 (cos 9' - cos 9'{0)) 5 (0' - <f>'(0)) (6) 

The joint probability of finding the rod with orientation u(0) in solid angle 
dQ'(0) at time and u(t) in dfl'(t) at time t can be written as 

p(fi'(t),t;fi'(0),0) 



2tt(1 - cos^o) 

E E e^[-i/n^ + i)i>|t|]yjr^(n / (o))y^(n'(t)) (7) 



oo oo 



n=l m=— oo 



The latter satisfies the normalization condition 

P (n'(t),t; n'(o),o)dn'(t)dn'(o) = 1 



cone cone 



where the angular integrals are taken only over the volume of the cone. 
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3 NMR framework for rotational diffusion 



At the beginning of this Sec. we recall the main facts from the general the- 
ory of spin-lattice relaxation by dipole-dipole interaction suggested by BPP 
and developed by Solomon [53]. The BPP-Solomon scheme is substantiated 
by more stringent Redfield's theory (see [17] for detailed presentation). This 
scheme is developed in the frame whose z axis is directed along the constant 
magnetic field. For the case of identical spins I the contribution to the spin- 
lattice relaxation rate constant due to rotational diffusion with the spectral 
density at a Larmor frequency ujl of the correlation function for spherical 
harmonics has the form (see VIII. 76 in [17] ) 

imxout = 374 " 2/ 2 (/+1) + j (2) (^l)} o) 

where 7 is the gyromagnetic ratio of the nucleus, / is their spin and h is the 
Planck constant. For non-identical spins / and S we have four equations (see 

7/ 2 7^ 2 S(S + l)x 



VIII.88 in @7J) 

(l/Tl 1 ) = 

\ ' 1 J rotat 



and 

( 1 /^)_=7?7ltf/(/ + Dx 

{-^'H-^ + l^H+^i)} < 10 > 



Here only two equations are written out explicitly because the other two can 
be obtained from them by mere changing of indexes |47) . To find the spectral 
densities J^(co), J^(co) and J^(co) we need to know the correlation functions 
G®(t) where i = 0,1,2 . 

As was stressed in the previous Sec. we set 8' and 0' to be polar angles defining 
the direction of the axis connecting protons and ty(8',<ft',t) = ^ (fi',i) to be 
the probability of the orientation of this axis in the direction Q' at time t. In 
the general case the axis of the cone can be tilted relative the magnetic field 
(z axis of the laboratory fixed frame) at an arbitrary angle ip. That is why 
at application to a realistic system the correlation function of internal motion 
in the cone < F^(0)F^(t) > internal has to be averaged over the orientations 
of cone axes with some overall distribution of the angles f(ip) characterizing 
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the system of interest. That is the correlation function G^(t) whose spectral 
densities are to be substituted into Q or ([ID]) has the form 

G«(f) =« F«(0)F«(t) > mt ernal>overall = 

- 1 # sin ^ /(V) < F»(0)F«(t) > internal (11) 
o 

To put it differently we assume that the distribution function /(?/>, A, u) (where 
ip,\,u are Euler angles for rotation of the dashed Cartesian frame x',y',z' 
relative the laboratory fixed one x, y, z) depends only on the Euler angle i/j, 
i.e., f(ip,\,uj) = f(ip)- It will be shown in Sec. 6 that only in this case the 
overall averaging provides the absence of cross-correlational functions with 
q ^ q', i.e., 

« F^(0)F«\t) > mt ernal>overall= &tfG<»\t) (12) 

where 5 nm is the Kronecker symbol. The latter requirement is crucial for the 
validity of the BPP-Solomon scheme [IT]. Further we consider four particular 

cases. 

a) . The cone axis is directed along the magnetic field for all cones in the system. 
It means that /(?/>) = S(ip) where S(x) is a Dirac 5— function (see Appendix 
B for technical details). This case is of little practical significance. However it 
provides direct comparison of the limiting case of our results with the textbook 
formulas from [I7J. Thus it serves as a test for the validity of the present 
approach from the theoretical side. Besides the formulas obtained in this case 
without superfluous complexities are further used in more involved cases as 
building blocks. That is why we denote the correlation functions G®(t) for 
this basic ones g®(t) = (t) f(^)=s(i/>) ■ This case is considered in 
details in Sec.4-Sec.5. 

b) . For the particular case of random isotropic distribution (unweighted aver- 
age) of cone axes relative the magnetic field we have f(ip) = 1. This case is 
considered in Sec. 6. 

c) . As an example of the case for the cone axes to be predominantly oriented 
along or opposite the magnetic field we consider the model function f(ip) = 
3 cos 2 ?/'. This case is considered in Sec. 6. 

d) . As an example of the case for the cone axes to be predominantly ori- 
ented transverse to the magnetic field we consider the model function f(ip) = 
3/2 sin 2 ?/;. This case is considered in Sec. 6. 
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From now and up to the end of Sec. 5 we consider the case a)., i.e., f(if>) = 
5(ijj) where 6(x) is a Dirac 5— function. In this case each cone axis is directed 
along the magnetic field and we need not distinguish the dashed Cartesian 
frame x', y', z' from the laboratory fixed one x, y, z. To retain the uniformity 
of designations for correlation function of wobbling in a cone we further use 
the dashed Cartesian frame and correspondingly the dashed spherical frame 
8', cf)' . We start from the general expression for the correlation functions of 
arbitrary order (see VIII. 13 in [17]) that in our case takes the form 

g ®(t) = G®(t) m=m =< F«(0)F«(t) > mternal = 

J J f« *(n'(t))F^(n'(o))p(n'(t),t;n'(o),o)dn'(t)dn'(o) (13) 

cone cone 



where i = 0,1,2 and the angular integrals are taken only over the volume of 
the cone. We need the correlation functions g^°'(t), g^(t) and g^ 2 \t) in order 
to substitute their spectral densities in Q or (flUi) . For our case of dipole- 
dipole interaction of two spins separated by the distance b they are defined by 
random functions F^°\ F^ and F^ [17] whose relationship with associated 
Legendre functions P^ (cos6) is known (see, e.g., Appendix c. in [55] ) 



„ ml 1 — 3 cos 2 9 2 „f 0i 

F I - 0) W = p = -p;Pf 



F« (Si) = ^'"^M = i/>» (cos 9 )exp W ) 

p(2) (n) = s..^exp(2^) = _l_ p < 2) (cos ff) exm4>) (14) 



We stress once more that in our particular case f(tp) = 5(tp) the (f2) in the 
laboratory fixed frame are identical to (Q 1 ) in the dashed (cone-related) 
frame to be substituted in ([TBI . 

Now we have to substitute ([TJ and (fbT|) into (fTBT) . We denote 

/i = cos 9' (15) 



so that fio = cos 9 and introduce the associated Legendre functions Pffl (fi) 

Y$> (n 7 ) = J-^eMirn^Pip (ji) (16) 
V l-ntin 



9 



which satisfy the orthogonality properties 



rti,ra 2 



(17) 



A'l) 



where 5„ iTO is the Kronecker symbol (S n>m — 1 if n — m and 5 n>m = other- 
wise). 

Making use of 1.12.1.12 and 1.12.1.9 from [56J (see Appendix A) respectively 
we obtain after straightforward calculations 



9 (l \t) 



1 + At 
6 6 



£exp[-^(^ + l)D |i 



n=l 



Tji 1 ) 
tin 



I i _ 9 , 2 { [« + 3)fig - l] P$ (aid) - ^oP^+i (Mo)}' (18) 



K + 3) 2 K-2)^ 



and 



(l-/i ) 2 (l + /i ) 



3 oo 



(i* + 3)*(itf -2) 



6 6 



n=l 
2 



£exp -z,> 2 + l)L>|t 



Pi 



(2) 



(19) 



The calculation of G^(t) requires the formula from [^8 



± 

K° n = Jdfi (3fi 2 - 1)P$ = 4z (l - 6z + Gz%) F u° n + 1; 2; z ) + 



/'0 



3^o (1 - 2z ) F (-!/«, z/° + 1; 3; z ) + 2z 2 F (-z/°, z/° + 1; 4; 



2f) 



(20) 



where 



2o 



1 - Mo 



(21) 



and P(a, 6 5 is a hypergeometric function. Making use of ( 120]) we obtain 



1 00 1 
^ E ex P + 1)P I t I] ^ 



(1 - Mo)& 6 fi 



(22) 
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We denote 



t. 



(m) 



v™{v™ + l)D 



(23) 



where m = 0,1,2. We obtain for the basic spectral densities 
and j( 2 \oS) of g^\t), g^ l \t) and g^\t) respectively 



J (0) M 



E 



-(0) 



— Ik°Y 

r(0) I «/ 



(24) 



and 



r> oo 

-V 

b 6 



-(1) 



n=l 1 + (uTn ) 



(1)\ 2 rr(l) 



X 



[1 + Mo) 



{ [(i/i + 3)jig - l] P$ (aid) - ^oP^+i (^o)} 2 (25) 



and 

.(I), , 2 ~ T< 2 ' (1 - gpf(l + ft;) 3 , ( 3) ,2 



These spectral densities enable us to calculate any of the spin-lattice relaxation 
rates ()9l-(IT0"l) for the particular case that the cone axis is directed along the 
magnetic field. 

We denote the correlation time r and the rotational correlation time r rota t 

87ia 3 r] Ana 3 !] r 

T = J^T Trotat = 3k^f = 6 ( ' 



4 Spin-lattice relaxation rate for nuclear pair with non-identical 
spins 



For the particular case that the cone axis is directed along the magnetic 
field we have for the spectral densities to be substituted in (19])- (1T0"]) are: 
J(°)(u;) = jW(cj), jW(w) ee jV>(u) and (to) ee j^(tu) where 
and j^(uj) are given by fT24|) - (!26|) . It is worthy to note that if we identify S 
with, e.g., 15 N from a nuclear pair of non-identical spins 15 N — H then we 
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actually need only the formula for (1/T/ 7 ) (see (fTUl) ) for the analysis of 

V / rotat 

experimental data. That is why further we restrict ourselves only with explicit 
writing out the formula for this quantity. The substitution of fl2^|) -( )26|) into 
(DUD yields 



(. 



b 6 



0\2 



■ x 



{ ((1-75/7/)^) 
3(l+/io) 



(1 + 75/7/) u{t 



1 + M + i)P 

{ + 3)tfj - l] P« (//o) - ^oP^ (/io)} 2 + 



-1 



3(l-^ ) 2 (l + /io) 3 



X 



2^ ) K + 3) 2 (z/ 2_ 2)2z/ 2 (z/ 2 + 



(2? 



This formula describes the spin-lattice relaxation rate from restricted rota- 
tional diffusion in a cone for the particular case of the cone axis to be directed 
along the magnetic field. The series in this formula is well convergent. That is 
why in practice it is sufficient to restrict oneself only by several initial terms 
in it. 



5 Spin-lattice relaxation rate for nuclear pair with identical spins 



For the case of identical spins the substitution of ( 1M1) - (1261) into (j9J) yields 



1 



r 7 4 n 2 /(J+l); \Tj rotat 
(1 + Ho) 



3 £kK + i)i 



1 



+ Kr) 2 \v\{y\ + 1)] 



^ 1} K + 3) 2 K-2)2 
1 



l + (2a; i r) 2 K(^ + l)]- 2 



x 



;i-/i ) 2 (i + /io) ; 



(29) 
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It will be shown later that in the limit of isotropic (6q — > tt) rotational diffusion 
this formula yields the well known result VIII. 105 from |47J. 



6 Arbitrary orientation of the cone axis relative the magnetic field 



In the general case of arbitrary tilted cone axis relative the magnetic field we 
need two frames (see Sec. 2 and Sec. 3). The laboratory fixed frame has the z 
axis directed along the magnetic field while the dashed (cone-related) frame 
has the z' axis directed along the cone axis. The angle between the cone axis 
and the magnetic field is ip. In (1131) we carry out internal averaging over the 
rotation in the cone in the dashed (cone-related) frame. That is why we need 
the transformation of the (f2) given by CHI) in the laboratory fixed frame 
into those in the dashed (cone-related) frame. 

As is well known a rotation of one frame relative the other is most conveniently 
described by Euler angles. We choose the angle ip as the first Euler angle 
(0 < ip < 7r). We denote two other Euler angles as A (that between the x-axis 
and the so-called iV-line (node-line) < A < 2tt) and u (that between the 
iV-line and the x'-axis < u < 2tt). The formula for transformation of the 
generalized spherical harmonics at transition from the frame {(p, 9} to that 
{0', 9'} obtained by rotation of the z axis by the angle ip is [58], [59] 

n 

P 2 (9) (cos#)exp(i#) = RtlW P2 S) (cos 9') exp [is^' + iq\ + isco] (30) 



where 



min(2— q,2— s) 
r=max(Q,—q—s) 



(2 + ?)!(2- a )! ( ^V +s+2r (. <A 4 " 9 - S - 2r ' 

1 cos — sin — (31) 



r\(2 - q - r)\(2 - s - r)\(q + s + r)\ \ 2) \ 2 

One can see that if the distribution function f(ip,X,u) characterizing a system 
of interest depends only on the angle ip, i.e., f(ip, A, uj) = f(ip) then at overall 
averaging 

7T 2-7T 2lT 

< •■• >overaii= ^3 J dip J d\ J du simjj f (ip, A, u) ... (32) 
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we have the factors 

2 IT 



J dX exp [i(q — q')X] = 2n5 qq i 



o 

2tt 



du exp [i(s — s')cj] = 2n5 s 



(33) 



(34) 



It is namely the identity ( 13"3"|) that provides the applicability of the formula 
( TT2|) . The latter is crucial for the validity of the BPP-Solomon scheme 
The overall averaging takes the form 



< ••• >overall= ~ J dtp SU1 1p f(lp). 




(35) 



and total (overall + internal) averaging is given by (11 ip in this case. The 
function f(ip) must be normalized so that 



dip smip f(ip) = 1 



(36) 



We denote 



h {2 



h {2 



*8 



dip sinipf(ip) (cos- j ( sin ^-) 



J dip smipf(ip) ( C0S ^J ( sin f) ( cos f 



o^> smipf(ip) 



cos- + sin- 



^■0 sin 3 ^ cos 2 ipf(ip) 



y sin ipf(ip) I - sin 4 r/> — cos 2 ^ sin 2 ^+ 



' 
sin — 
, 2 . 



COS 2 



, si 4) + ( co 4) 
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/i (0) 



sin 



cos • 



1 > 2 

c£0 sin ipfiip) ( cos2 ^ — 2 s ^ n2 ^ 



/i (0) - h (1) 



h < § ) = /# sin 5 ^) 



(37) 



After lengthy but straightforward calculations we obtain 



(38) 
(39) 
(40) 



where j(°\u), and j^(oj) are given by (1!Mj) - f[2"6T) . In designations (1271) 

the explicit form of the basic spectral densities j(°\u)), and ^ 2 '(w) is 



J (0) M 



2r 



E 



& 6 (l-^o)^ 1 M°Vo( z/ o + 1 ) 



K} 1 



1 + 



[COT 



\2 I- 1 



<k + i)n 



(41) 



(1) _ 2r(l + ^) ~ K + 3)- 2 (i/jj-2)- 2 w 
' (W) - _ S^fe ^, + 1) 

{[(z/J + 3)^ - l] P« (/io) - i/^ioP« +1 (// )} 5 



, (2) 2r(l- / i ) 2 (l + /,o) 3 ^ (^ + 3)- 2 (i^-2) w 



2 



(42) 



b 6 



(3) 



1 + 



n=l 
2 n-1 



M 2 V 2 ^ 2 + 1) 



+ 1)]' 



(43) 



For practical application of the theory we consider three cases. 
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1. For the particular case of random isotropic distribution (unweighted aver- 
age) of cone axes relative the magnetic field we have f{ip) = 1. In this case 

we obtain: hg) = V 15 ! ^(i) = V 5 ! h U = 4 A kjjj = 4/15; ^(i) = 8 / 15 ! 
h§ = 4/5; fcgj = 2/5; Jig} = 4/15; fcgj = 16/15. 

2. As an example of the case for the cone axes to be predominantly ori- 
ented along or opposite the magnetic field we consider the model function 
f(ip) = 3 cos 2 ?/'. In this case we obtain: h$. = 1/35; h^l = 1/7; /tSj = 44/35; 

h$ = 12/35; h$ = 28/35; hfy = 20/35; = 22/35; h$ = 12/35; 

= 16/35. 

3. As an example of the case for the cone axes to be predominantly ori- 
ented transverse to the magnetic field we consider the model function f(ip) = 
3/2sin 2 ^. In this case we obtain: ftiS = 3/35; h?X = 8/35; h$) — 20/35; 

= 8/35; h$ = 14/35; hfy = 32/35; h$ = 10/35; h$ = 8/35; = 
48/35. 

Making use of fj4Tl) - fj43l and explicit values of h$ for these three cases enables 
us to calculate the spectral densities fl38l - fj40l to be inserted in BPP-Solomon 
formulas fl9l)- f|T0l . We do not write out explicitly the corresponding expressions 
for the spin-lattice relaxation rate to save room. However in the next Sec. we 
plot the corresponding dependencies of spin-lattice relaxation rate for nuclear 
pair with non-identical spins for all three cases. 



7 Results and discussion 

For 15 N — H nuclear pair of non-identical spins, we identify S with 15 N 
nucleus and / with H one. Thus 75 = —2712 rad s' 1 Gauss' 1 and 77- = 
26753 rad s' 1 Gauss' 1 so that 75/7/ = —0.101372. For 13 C — H nuclear pair 
75 = 6728 rad s _1 Gauss' 1 so that 75/7/ = 0.251486. The righthand side in 
the formula (128]) depends on the parameters characterizing the nuclear pair 
(namely on the gyromagnetic ratios 75 and 7/ of our pair of non-identical 
spins). To plot the spin- lattice relaxation rate with the help of (T28]) one has 
to choose the particular nuclear pair explicitly. That is why to be specific we 
choose the 15 N — H nuclear pair of non-identical spins. It should be mentioned 
that both (1281) and (1291) have uncertainties at: 6$ = tt/A because v\ = 2.0000; 
#0 = vr/2 because v\ = 2.0000; 6 = 3vr/4 because v\ = 2.0000; 6 = 175° 
because v\ = 2.0000. However these uncertainties are isolated, and can be 
safely ignored. 
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In Fig. 1 the spin-lattice relaxation rate for identical spins obtained with the 
help of ( )29|) for the case of cone axes directed along the magnetic field is 
depicted as a function of the cone half- width 8q at different values of Larmor 
frequency. The value ult = 0.1 for the upper curve is within the range of 
validity of the extreme narrowing limit lolt « 1. From this Fig. one can 
see that for the case of isotropic (8q = n) rotational diffusion in the limit of 
extreme narrowing the corresponding curve tends to the value 0.333.... Thus 
the formula fl29|) yields 



lim 



7T \TJ i h 2 I(I+ 1 




rotat 



a 



1 

3 



(44) 



that taking into account ff27|) coincides with the well known formula VIII. 106 
from 



1) = *g I{I + 1) g* 

T\ J rotat 3fc B T 



In Fig.2 the spin-lattice relaxation rate for non-identical spin pair 15 N — H 
obtained with the help of lj3"51)-(l4"01) for the case of random isotropic distri- 
bution (unweighted average f(ip) = 1 taking place in powders) of the cone 
axes relative the external magnetic field is depicted as a function of the cone 
half-width 8q at different values of Larmor frequency. The reduced curves 
[b 6 / (T-fhs h2s ( s + !))] (l/^i 1 ) for 13c - H nuclear pair are found to be 
very similar to those for 15 N — H one. For this reason data for 13 C — H nuclear 
pair are not presented to save room. 

In Fig. 3 the dependence of the spin-lattice relaxation rate on the cone half- 
width #o for the case of non-identical spin pair 15 N — H is depicted for different 
distributions of the cone axis relative the external magnetic field: for the case 
of cone axes directed along the magnetic field (f(ip) = S(ip)); for the case of 
the cone axes to be predominantly oriented along or opposite the magnetic 
field (/(VO = 3cos 2 ^ that may be relevant for, e.g., liquid crystals); for the 
case of random isotropic distribution (unweighted average f{ip) = 1 taking 
place in powders); for the case of the cone axes to be predominantly oriented 
transverse to the magnetic field {f{ip) = 3/2 sin 2 ijj that may be relevant for, 
e.g., liquid crystals). 

As is well known in liquids (for which isotropic rotation is relevant) T\ ordi- 
narily decreases with increasing viscosity, in some cases reaching a minimum 
value after which it increases with further increase in viscosity [J7]. The vari- 
ation of viscosity is caused by temperature T. For ordinary isotropic rotation 
correlation times are functions of temperature r n = r (tj(T)) [n(n + 1)] _1 (see 



17 



VIII. 97 in [17] ) and the spin-lattice relaxation rate (l/Ti) rotat (r(T)) can have 
a maximum as a result. In the present model we have one more option. Sub- 
stitution of (JZTJ into ([23]) yields 



(m) 



r(v(T)) 



(45) 



n 



where m — 0, 1, 2. Thus correlation times are functions of temperature and of 
the cone half-width #0 (via the values v™). The results obtained testify that 
variation of each of these parameters (T or # ) can produce a maximum in the 
corresponding dependence of the spin-lattice relaxation rate. The dependence 
in Fig.l exhibits a maximum at 6 w 100°. Fig. 2 and Fig. 3 show similar 
maximums for non-identical spin pair at 6q 7r/2. These maximums are 
similar to those at #0 ~ 7r /2 given by the model-free approach. 

Indeed let us consider the most typical for practice case of the model-free 
approach when the overall motion of a macromolecule is considerably slower 
than the internal motion. In this case the expression of the model-free approach 
for the relationship of the spin-lattice relaxation rate with the order parameter 
S 2 and effective correlation time r e is given by equation (37) from [13] 



where a and b are constants independent on spatial configuration accessible for 
internuclear vector. For wobbling in a cone the relationship between the order 
parameter of the model-free approach and the cone half-width 6 is given by 
equation (A3) from [13] 



while that for the effective correlation time is given by equation (A4) from 




(46) 



S = - cos 6*o(l + cos 6*o) 



(47) 



[43j 



1 



cos 6 (1 + cos 6 ) 2 {In [(1 + cos o )/2] + 



D w (1 - S 2 ) 



(1 



cos 9 )/2}/ [2 (cos 9 - 1)] + (1 - cos# )(6 + 8 cosfl, 




(48) 



18 



where D w is the diffusion coefficient. We denote 



Then we obtain the dependence of the spin-lattice relaxation rate on the cone 
half-width 6 Q 



-— = - cos 2 O (1 + cos 9 ) 2 + c\ cos 2 6 (1 + cos 6 ) 2 {In [(1 + cos 6 )/2] + 



The dependence of this spin-lattice relaxation rate on the cone half-width 6 
is depicted in Fig. 6 at several values of the parameter c. One can see close 
qualitative similarity of the curves with our results. 

The BPP-Solomon scheme requires essential extent the isotropy (but not to- 
tal isotropy) for its validity. In Sec. 6 we show that in application to rotational 
diffusion in a cone it remains valid for systems with a distribution of cone 
axes depending only on the tilt relative the magnetic field but otherwise being 
isotropic. This residual isotropy provides the requirement ( 155|) that is neces- 
sary for the absence of cross-correlational functions (IT2|) . The latter in turn 
is crucial for the validity of of the BPP-Solomon scheme. We consider the 
aforesaid a so important issue that would like to reiterate it in other words 
with complete definiteness. We develop the theory for the general case of arbi- 
trary orientation of the cone axis relative the magnetic field (laboratory fixed 
frame). We show that when the cone axis is tilted at an arbitrary angle ip 
to the magnetic field but otherwise is oriented isotropically then at overall 
averaging the crucial requirement for the validity of the Bloemberger, Pur- 
cell, Pound - Solomon scheme ( !T2|) (that of the absence of cross-correlational 
functions with q ^ q') is retained. Thus we explicitly prove the consistency 
of combination of the textbook formulas for the BPP-Solomon scheme with 
rotational motion in a cone. 

Practical applications of our results depend on the choice of a distribution 
function for overall averaging over all orientations of the cone axis with respect 
to the laboratory fixed frame (over the angle ip). This distribution function 
f{ip) is a characteristic of the system of interest. Under the assumption of 
isotropic random orientation of cone axes relative the laboratory fixed frame 
the results obtained can be applied to powders. We consider this case of un- 
weighted average over all orientations of the cone axis with respect to the 



(1 -cos#o)/2}/[2(cos# - 1)] + (1 - cos 0„) (6 + 8 cos 6 t 




(50) 
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laboratory fixed frame (f{ip) = 1) in Sec. 6 and provide corresponding data 
for the spin-lattice relaxation rate in Fig.3. Also in Sec. 6 we consider a model 
example of predominant orientation of the cone axis along or opposite the 
magnetic field (f(ip) = 3cos 2 ip) and that of their predominant orientation 
transverse to the magnetic field (f(ip) = 3/2 sin 2 ip). The results for these 
cases may be relevant for, e.g., liquid crystals. For both of these cases we also 
provide corresponding data for the spin-lattice relaxation rate in Fig.3. The 
particular case of the cone axis directed along the magnetic field {f{ip) = S(ip) 
where S(x) is a Dirac delta-function) is of little practical significance. However 
this case provides direct comparison of the limiting case of our formulas with 
the textbook formulas from [37]. Thus it serves as a test for the validity of our 
approach from the theoretical side. 

Fig.3 shows that the results for practically relevant cases of specific distribu- 
tions of cone axes relative the magnetic field are qualitatively similar to those 
for the model case of cone axes directed along the magnetic field. However 
there is some quantitative difference. From Fig.3 one can see that 

\Ti) f^ =S ^ \Ti) /(^) = 3 cos 2^ (ji) f(f)=l (ti) /(</,)=3/2sin 2 V 



Thus the more is the contribution of the orientations of the cone axes in the 
system transverse to the magnetic field (i.e., the more is the probability of 
such orientations) the more efficiently spin-lattice relaxation proceeds. This 
dependence is weak but well-defined. To prove it explicitly we plot in Fig.4 
the the dependence of the spin-lattice relaxation rate on the tilt angle of 
the cone axis relative the magnetic field for the cone with the half-width 
#o = 55° as an example. This case corresponds to the distribution function 
f(ip) = 5(ip — ip ) [cos (-^o/2)] 1 (see Appendix B for technical details). One 
can see that when the cone axis is transverse the magnetic field the spin-lattice 
relaxation proceeds most efficiently In Fig. 5 such dependence is plotted for 
different values of the cone half- width 6q. 

The rigorous quantum-mechanical treatment presented in the main text of the 
paper does not allow simple physical interpretation of the results obtained. To 
attain the latter objective we invoke to the results of the classical analogy pre- 
sented in Appendix C. Such simplified approach does not yield precise results 
but enables us to get more clear physical perception of the phenomenon under 
consideration. We conclude that the essence of the maximum in the depen- 
dence of the spin-lattice relaxation rate on the cone half-width 8q originates 
from two factors: 

1. the angular dependence of the spin-lattice relaxation rate 
— oc cos 2 8 sin 2 9 
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and 

2. the averaging of the spin-lattice relaxation rate over the rotational diffusion 
in the cone with half-width O 



Combination of these two factors leads to the required non-monotonic behavior 
(see Fig. 7). Regretfully further simplification of the physical picture underlying 
the phenomenon seems hardly possible. It is hampered namely by the fact that 
the phenomenon can not be attributed to a single reason but originates from 
the combination of two of them. 

Thus the model yields an unexpected result: too much spatial freedom is 
not good for spin-lattice relaxation rate. This result is obtained from rigor- 
ous quantum- mechanical treatment, is loosely supported by its rough classical 
analogy and besides is well corroborated by the model-free approach. The re- 
striction of free isotropic motion (9 = ir) to smaller cone half-widths leads 
to the increase in spin-lattice relaxation rate with a maximum in the region 
of 6q mx 90° -T- 100°. Only further restriction of the cone semi- angle to the 
values 6q < 0™ ax leads at last to efficient decrease of the spin-lattice relax- 
ation rate. We can reformulate the above mentioned physical interpretation of 
this phenomenon in other words. The maximum arises at 6$ ~ 6™ ax because 
in this case the internuclear vector has maximal probability to be transverse 
the external magnetic field. As one can see from Fig.4 and Fig. 5 such orien- 
tation maximizes the spin-lattice relaxation rate. Increase of 6q from 6™ ax to 
the limit of free isotropic motion (8q = 7r) decreases the probability for in- 
ternuclear vector to be transverse the external magnetic field and thus leads 
to the decrease of the spin lattice relaxation rate. Taking into account that 
wobbling in a cone is the exactly tractable model we can not attribute this 
result to any approximations. Also the model captures the effect of restriction 
for the motion in an adequate and non-trivial manner. Thus the phenomenon 
can not be model dependent. We conclude that the phenomenon should be 
of general character, i.e., any correct model of restricted rotational diffusion 
should exhibit a non-monotonic dependence of the spin-lattice relaxation rate 
on the accessible space volume for the internuclear vector. The maximum for 
this dependence should take place for those configurations for which the prob- 
ability of the inter-nuclear vector to be transverse the external magnetic field 
is maximal. 

This result can manifest itself in experiments where addition of some substance 
creates efficient steric hindrance to rotating internuclear vector of the reporting 
nuclear spin pair. Its accessible space volume initially (without this substance) 
ought to be large enough so that one can approximate it by the limit of 
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free isotropic motion 9q — > n. Also one should take special concern that the 
additional substance does not alter the micro-viscosity for the environment 
of the reporting spin pair appreciably. It can be achieved if the size of the 
molecules for the added substance is not small but rather commensurable 
with those containing the reporting nuclear spin pair. Only in this case one 
can obtain the required condition that variation of the concentration C of 
the added substance touches upon primarily v™(C) via the accessible space 
(cone half-width) while micro- viscosity remains mainly constant rj(C) ~ const. 
Otherwise we are in danger that variation of the correlation times 

r M = r( V (C)) 
n *C(C)(^(C) + 1) 



is attributed namely to 77(C) rather than to u™(C). In the latter case one actu- 
ally probes the dependence 1/Ti vs. micro- viscosity 77(C) with non-monotonic 
behavior of quite the same nature as that for the usual dependence 1/Ti vs. 
temperature T stipulated by the dependence n(T). Provided all precautions 
are taken into account and the above condition is satisfied one may hope that 
the dependence 1/Ti vs. cone half-width 8q can be probed in the experiment. 
Then basing on the results obtained in the present paper (Fig.l, Fig. 2 and 
Fig. 3) one can expect that the increase of the concentration of the additional 
substance should lead to the increase of the spin-lattice relaxation rate up to 
a maximum with further decrease. The results of our rigorous treatment will 
be more accurate for the description of such experiments than those of the ap- 
proximate approach jl3], [51]. Indeed Fig. 6 shows that the latter predicts the 
maximum strictly at one value of the cone half- width, namely at 9™ ax « 85°. 
On the other hand Fig.l, Fig.2 and Fig.3 show that within the framework of 
the rigorous treatment the position of the maximum is a variable value from 
the range 80° < 6™ ax < 180°. It depends on the type of the nuclear spin pair 
(homonuclear or heteronuclear) , on the Larmor frequency or more exactly on 
the value of ojlt, etc. The dependence of the maximum position on the type 
of the labeled atom ( 15 iV or 13 C) is found to be very small. Thus we antic- 
ipate that the results of the present approach leave much more freedom for 
describing experimental data. 

We conclude that wobbling in a cone is the exactly tractable (i.e., amenable to 
full analytic treatment within the range of validity of BPP-Solomon scheme) 
model for description of the nuclear spin-lattice relaxation rate. This fact 
makes it a unique one among all other models for restricted rotational diffu- 
sion. 
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8 Appendix A 



Here we present two known mathematical formulas 1.12.1.12 and 1.12.1.9 for 
the the associated Legendre function Pj^ (x) from the table of integrals 
The formula 1.12.1.12 is 



J dx x(l - x^^Pj? 



:i - x 2 )^ 2 



(i/±/i + 2)(i/t/i- 1) 
{ \{v ± fi + 2)x 2 - l] P>? ( x ) + (ji-u- l)xP» +l (x) } 



The formula 1.12.1.9 is 

dx{l- x 2 yl 2 P» (x) = 1 X ) P^ (x) 

J (U — fJL){U + fJL + 1) 



9 Appendix B 



The polar angle ip at operations with the distribution function f{ip) imposes 
some peculiarities in treating the case of Dirac 5— function 

fW = f ) l0 M - V>o) 

cos ('00/2) 



We stress that the case of the cone angle oriented along the magnetic field 
f(ip) = S(ip) considered in Sec. 4- Sec. 6 is a particular case of this distribution 
function corresponding to the value ipo = 0. First let us prove the normaliza- 
tion requirement (I36I) . We have (taking into account that 5(2z) = S(z)/2) 

~ J dip sinip 5(ip-ip ) = 2 J d ( — J sin ^ cos ^5 \ ' " 

o o ^ ' \ 



2 



tt/2 tt/2 



/ dx sin x cos x 5 (x — y j = — j d(cosx) cosx 5 (x — y j 

^ ' ^ ' 

J dy y 5 ^arccosy - y^J = J dy y 5 (^y - cos y^J 



00 

cos — 
2 



This calculation serves as a model for operations at calculating the values of 

An) 



hfa\ in (|37|) with the distribution function /(0) = 5(ip — O ) [ cos (0o/2)] 1 - 
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The general rule takes the form 



1 



TT 



cos (ip /2) 



dtp sinip S(ip- 4> )q(?P) = 2g(^ ) 



o 



where q{ip) is an arbitrary function of if). 

At %I)q = we obtain = 2, = 2 and h^l = 2 while h^A = at 
m 7^ n. Substitution of these values into (I38l) - (|40l) yields J^ 2 \u) = j^^u), 
J(i)( w ) = jW(co) and J(°)(w) = j (0) (w) as it must be. 



10 Appendix C 



The rigorous quantum-mechanical treatment described in the text of the 
present paper does not enable us to gain insight in physical interpretation 
of the results obtained. To attain the latter objective we present here its clas- 
sical analogy that is much more simple and not at all precise but nevertheless 
rather illustrative. Let us consider the simplest case of two identical spins (nu- 
clei) with magnetic moments fx which we assume for simplicity to be oriented 
along the external magnetic field Hq directed along the axis z of the spherical 
frame so that fx = fie z . The internuclear vector (with the absolute value r) is 
tilted by the angle 9 relative the external magnetic field. We denote n the unit 
vector for internuclear one. Then the magnetic field produced by one spin in 
the location of the other is 

h oc -g [jt — (jx ■ n) n] 



The spin-lattice relaxation rate for classical treatment of dipole-dipole inter- 
action is given by the following expression (see, e.g., Appendix E in [60J) 

1 °° 

— = 1 2 J dr cos (u l t) [< h x (0)h x (T) > + < h y (0)h y (T) >] 



where 7 is the gyromagnetic ration and ujl = 7 Hq is the Larmor frequency. 

For our diffusion in a cone model we have to average the latter value over the 
rotational motion in the cone with half- width 8 . Thus we want to calculate 
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the value 




1 ' cone 



OC 



1 — cos 9, 



1 



'o 



( dO s\n6 — 



The unit vector n in the spherical frame is 
n = sin 9 cos (pe x + sin 9 sin (fie y + cos 9e z 

so that 

(/J • n) — /i (e 2 • n) — /i cos 9 

One can see that both h x oc cos 9 sin 9 and /i y oc cos 9 sin so that 
J- oc cos 2 sin 2 9 



The origin of the latter relationship can be also trivialized as follows. The 
energy for the dipole-dipole interaction of two magnetic moments is 

U oc 3 cos 2 9 - 1 

Then the force of interaction (this notion is quite viable in our classical analogy 
of the quantum-mechanical treatment) is 

F oc — — oc cos 9 sin 9 
d9 

The spin-lattice relaxation rate is determined by the correlation function of 
the force that yields the above relationship. Thus we obtain that the angular 
dependence of the spin-lattice relaxation rate on the cone half-width 9 is 



i 




'0 







Integration yields the final result 




cos' 




This dependence is depicted in Fig. 6. 
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Fig. 1. Spin-lattice relaxation rate for nuclear pair of identical spins from rotational 
diffusion for the case of cone axes directed along the magnetic field (eq. (f29|) ) as 
the function of the cone half-width 8q (in degrees) at different values of the Larmor 
frequency: ult = 10 (thin dots), w^r = ujlt = 1, lolt = 0.1 (thick dots). 

For the case of free isotropic rotational diffusion (cone half-width 9q — > ir) the 
textbook value 0.333... ( i.e., formula VIII. 106 from |47j ) is obtained as it must be 
(see (eq. flU]) and text below). 
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Fig. 2. Spin-lattice relaxation rate for 15 7V — H nuclear pair of non-identical spins 
from rotational diffusion for the case of random isotropic distribution (unweighted 
average f(ip) = 1 taking place in powders) of the of the cone axis relative the labo- 
ratory fixed frame (external magnetic field) as the function of the cone half-width 
#o (in degrees) at different values of the Larmor frequency: lolt = 10 (thin dots), 
ult = 10 - 5 , ujlt = 1, ujlt = 0.1 (thick dots). 



31 



(1/TiW < 7(7?7?* 2 tS(S+1)) 
0.4 r 



x* - ■: i 



0.3 



0.2 



0.1 



f(iA)=<5(iA) 



f(iA)=3 cos 2 iA 
f(«=l 



50 



100 



150 



^degrees) 



Fig. 3. Spin-lattice relaxation rate for ./V — H nuclear pair of non-identical spins 
from rotational diffusion as the function of the cone half-width 9q ( m degrees) at 
lolt = 0.1 for different distributions of the of the cone axis relative the laboratory 
fixed frame (external magnetic field): for the case of cone axes directed along the 
magnetic field (f(ip) = for the case of the cone axes to be predominantly 

oriented along or opposite the magnetic field (f(ip) = 3 cos 2 ip that may be relevant 
for, e.g., liquid crystals); for the case of random isotropic distribution (unweighted 
average /(VO = 1 taking place in powders); for the case of the cone axes to be 
predominantly oriented transverse to the magnetic field (f(ip) = 3/2 sin 2 ^ that 
may be relevant for, e.g., liquid crystals). 
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Fig. 4. Spin-lattice relaxation rate for N — H nuclear pair of non- identical spins 
from rotational diffusion in the cone with the half-width 9q = 55° as the function 
of the tilt angle ipo of the cone axis relative the magnetic field (corresponding to 
the distribution function f(ifi) = 5(ip — V>o) [cos (V>o/2)] _1 ) for different values of the 
Larmor frequency: w^r = 10~ ' 5 (thick line); ujlt = 10 ' 2 ; w^r = 10 ' 9 ; ujlt = 10 
(thin line). 
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Fig. 5. Spin-lattice relaxation rate for N — H nuclear pair of non- identical spins 
from rotational diffusion in the cone at ujlt = 0.1 as the function of the tilt angle 
V'o of the cone axis relative the magnetic field (corresponding to the distribution 
function f(ip) = 5{ip — ipo) [cos (tpo/2)]~ 1 ) for different values of the cone half-width: 
6 = 35° (thin line); 6 = 55°; 9 Q = 95°; 6 = 130°; 9 Q = 170° (thick line). 
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Fig. 6. Spin-lattice relaxation rate for the model- free approach (eq. (|50p ) as the 
function of the cone half- width 9q f° r different values of the dimensionless parameter 
c (eq. gSD): c = 50 (thick line); c = 70 ; c = 90 ; c = 110 (thin line). 
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Fig. 7. Angular dependence of the spin-lattice relaxation rate on the cone half- width 
6q for classical analogy of the rigorous treatment of diffusion in a cone model (see 
Appendix C). 
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